A Monte Carlo Calculation of the Pionium 
Break-up Probability with Different Sets of Pionium 
Target Cross Sections 

C Santamarina^, M Schumann', L G Afanasyev§ and T Heim' 

f Institut fur Physik, Universitat Basel, 4056 Basel, Switzerland 

| Institut firr Theoretische Physik, Universitat Basel, 4056 Basel, Switzerland 

§ Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 

Abstract. Chiral Perturbation Theory predicts the lifetime of pionium, a hydrogen- 
like tt + tt~ atom, to better than 3% precision. The goal of the DIRAC experiment at 
CERN is to obtain and check this value experimentally by measuring the break-up 
probability of pionium in a target. In order to accurately measure the lifetime one 
needs to know the relationship between the break-up probability and lifetime to a 
1% accuracy. We have obtained this dependence by modeling the evolution of pionic 
atoms in the target using Monte Carlo methods. The model relies on the computation 
of the pionium-target atom interaction cross sections. Three different sets of pionium- 
target cross sections with varying degrees of complexity were used: from the simplest 
first order Born approximation involving only the electrostatic interaction to a more 
advanced approach taking into account multi-photon exchanges and relativistic effects. 
We conclude that in order to obtain the pionium lifetime to 1% accuracy from the 
break-up probability, the pionium-target cross sections must be known with the same 
accuracy for the low excited bound states of the pionic atom. This result has been 
achieved, for low Z targets, with the two most precise cross section sets. For large Z 
targets only the set accounting for multiphoton exchange satisfies the condition. 
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1. Introduction 

Pionium is the hydrogen-like electromagnetic bound state of a 7r + 7r~ pair. Its lifetime 
is determined by the strong 7r + 7r~ — > 7r°7r° annihilation in the low relative momentum 
regime where Chiral Perturbation Theory applies. The theory predicts the pionium 
lifetime value of (2.9 ± 0.1) x 10~ 15 s [1J. Experimentally, this value is being measured 
in a model-independent way by the DIRAC experiment at CERN [2]. 

Under the experimental conditions of DIRAC the pionic atoms are created in the 
inelastic scattering between 24 GeV/c protons, supplied by the PS at CERN, and the 
nuclei of the target [Hill], a chemically pure material of well-determined thickness. 

A pionic atom propagating at relativistic speed inside the target interacts mainly 
electromagnetically with the target atoms. The electromagnetic cross section is on the 
order of a megabarn. These interactions can lead either to a transition between two 
bound states of pionium or to a dissociation (break-up). The scattering with the target 
atoms competes with the process of annihilation since a shorter pionium lifetime leads 
to a smaller break-up probability. 

The main thrust of DIRAC's experimental technique is the detection of 7r + 7r~ pairs 
with very low relative momentum. From the low-momentum part of the spectrum, the 
break-up probability of pionium is then determined. 

The goal of this study is to establish the theoretical dependence of the break-up 
probability on the lifetime with a very high accuracy. We modeled the dynamics of 
pionium in different targets taking into account a large number of atomic shells that 
become populated as the pionic atom evolves in the target. 

We have performed the calculations by using three sets of pionium target interaction 
cross sections. The original studies by Afanasyev and Tarasov [H] made use of the Born 
approximation and pure electrostatic interaction. In our work we have applied new 
corrected cross sections that take into account relativistic effects |H] and the multiphoton 
exchanges [Zj. In addition, these cross section sets include a more accurate description 
of the target atomic form factor. As an error in the break-up probability translates 
directly to an error in lifetime, we have checked whether the results obtained in the new 
calculations significantly deviate from the previous ones. We conclude that the magnetic 
and relativistic corrections together with the target form factor choice are non-negligible 
mainly for the small Z targets, whereas the multiphoton exchange should be considered 
for the large Z ones. 

2. Monte Carlo Simulation of Pionium in the Target 

As we have noted in the introduction, to determine the pionium lifetime, the experiment 
DIRAC needs a precise theoretical calculation of the break-up probability P^ of pionium 
due to the electromagnetic interaction with target atoms. This calculation can be done 
by means of a Monte Carlo transport code that simulates the evolution of pionium 
from its creation to either its annihilation or its break-up under the given experimental 
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conditions. 



2.1. Pionium Production 

The pionic atoms are formed as a consequence of the Coulomb final state interaction 
of two oppositely charged pions. These pions are created in an inelastic collision of 
a 24 GeV/c proton and one of the target nuclei. The pionic atom is described by six 
quantities, accounting for the six degrees of freedom of a two body system. A particularly 
convenient choice of these quantities is given by the laboratory momentum of the center 
of mass of the atom P and the quantum numbers of the created bound state n, I, and 
m in the spherical coordinate representation \. 

The probability of pionium being created is given by jS| 



- (Ztt) |VW(<JJ| 



(1) 

P=q=P/2 



dP M dpdq 

The two terms on the right-hand side of the equation illustrate the final state interaction 
mechanism. The rightmost factor is the doubly inclusive cross section of 7r + and 7r~ pairs 
at equal momenta (p = q) without considering the final state interaction, as indicated by 
the superscript 0. The subscript s means that only pions created from direct hadronic 
processes and decays of resonances with a very short lifetime are considered, because 
the Coulomb interaction of pions from long-lived sources (e.g. T], Kg and A) is negligible 
and hence they do not contribute to pionic atoms production. The effect of the final 
state Coulomb interaction is to create a bound state with quantum numbers n, I, and 
m; it is given by the squared wave function at the origin. 

The doubly inclusive cross section can be obtained from the direct measurements 
of time correlated 7r + 7r _ pairs in DIRAC, according to the following reasoning: 

• The final state Coulomb interaction for short-lived sources is given, as in the case of 
the creation of a bound state, by a multiplicative factor depending only on Q, the 
magnitude of the relative momentum between the two pions. This is the so-called 
Coulomb or Gamow factor 9\ 

Aa ' -MQ& MQ) - , ™dL . P) 



dpdq uv ^dpdg' UK ^ J 1 - e ~ 27rM WQ 

where a is the fine structure constant. 

The contribution to the doubly inclusive cross section of pairs containing at least 
one pion from a long-lived source, ui(P), can be calculated with a hadron physics 
Monte Carlo simulation. In our case we have used FRITIOF6 JU]- This function 
has been shown to depend only on P |TT] . the magnitude of the total momentum 
of the pion pair. Taking this into account together with (J2J we find 

do" dcr, do~i . dcr° da 



dpdq dpdq dpdq dpdq dpdq 

thus relating a and of. 
\ Also parabolic quantum numbers have been used elsewhere jS] 
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Figure 1. The left graph shows the laboratory momentum magnitude distribution 
and the right graph the angular distribution with respect to the proton beam axis for 
low relative momentum tt + it~ pairs. 



Finally, we have found that the P-dependence of the doubly inclusive cross section 
is not correlated to Q, given that Q -C 30 MeV/c. 

These findings allow us to relate the P-dependence of a and a® by 

d ,Q~2McV/c d ^ 

K Jo (1 "^ (P)) d^ dQ ' (4) 



dpdq 



P=q=P/2 



where the P distribution is obtained from the direct measurement of the laboratory 
momentum of low relative momentum tc + tt~ pairs in DIRAC. In figure ^ we show the 
distribution of the magnitude of the momentum P and the angular distribution relative 
to the proton beam axis for low relative momentum 7r + 7r~ pairs. 

The initial quantum number distribution depends on the value of the wave function 
at the origin. It has been shown 12^ that the effect of the strong interaction between 
the two pions of the atom significantly modifies |V>nZm(0)| i n comparison to the pure 
Coulomb wave function. However, the ratio between the production rate in different 
states has been demonstrated to be kept as for the Coulomb wave functions ^3]. Then, 
considering that the Coulomb functions obey 



if I ^ 0, 




cM, 12? (5) 



we note that only S states are created, according to a 1/n 3 distribution. 

Another quantity to be specified is the position where the proton-target interaction 
took place (R). This is also the position where the atom was created. Since the 
target thickness is chosen much smaller than the nuclear interaction length of the target 
material, the tt + tt~ atoms are supposed to be uniformly generated all across the target 
thickness. The position in the transverse coordinates is unimportant, but it can be 
generated according to the beam profile characteristics, too. 
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We have now gathered all the information needed to simulate the creation of 
pionium atoms with a center of mass momentum P according to the experimentally 
measured momentum and angular distributions, at a position R uniformly distributed 
through the target, and in an initial S'-wave state with the principal quantum number 
n distributed according to a 1/n 3 distribution. 



2.2. Pionium Annihilation 

Once an atom has been created in its initial state, specified by P, n, I, and m, its 
dynamics are those of a free system that can either be annihilated, mainly via the 
7r + 7r~ — > 7r°7r° channel, or be electromagnetically scattered by one of the target atoms. 

The strong interaction decay to two neutral pions determines the lifetime and is 
related to the Oq — scattering lengths difference and to the wave function at the origin 
by dU 



u ^ -a^(l + S v ) O) , (6) 

where M n and M n o are the masses of the charged and the neutral pion, respectively, and 
Sr is the correction to next-to-leading order (5r = 0.058) that includes the effect of the 
strong interaction between the two pions. Using Chiral Perturbation Theory, Colangelo 
et al PP have been able to calculate the most precise value of the scattering lengths 
difference to date (a° — a% = 0.265 ± 0.004). Employing this value in © yields 

t 100 = (2.9±0.1)-10- 15 s. (7) 

Note, however, that due to (0) pionium may only decay from S states. Moreover, the 
lifetime of any S state is related to the lifetime of the ground state, by 

r„oo = n 3 r. (8) 

For the purpose of simulating pionium in the target, we shall from now on refer to r as 
the pionium lifetime. 

Hence, the probability for a tt + tt~ atom to annihilate per unit length, after the 
Lorentz boost transformation to the laboratory system, is given by 

1 ( J^Zl -f/-n 

PnL = W ,= \ PnH Xt ;- U ' (9) 
A mm I other cases, 

where is the annihilation mean free path. 
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2. 3. Electromagnetic Interaction of Pionium with the Target 

The electromagnetic pionium-target interaction of a pionic atom in an initial nlm bound 
state can induce a transition to another n'l'm' bound state. The probability of such an 
interaction per unit length is given by 

n'l'm' _ n'l'm 1 n\ 

Pnlm a nlm \ / 
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where p is the target density, A its atomic weight, No is the Avogadro number, and 
Onim 1 ' are the discrete (bound-bound) transition cross sections. 

The break-up mechanism is analogous to the discrete one; the break-up probability 
per unit length of an atomic bound state nlm is given by 



Pnlm \br A ®nlm (11) 

A nlm A 

where a^ lm is the break-up (ionization) cross section. 

Finally, the total cross section gives the probability of an atom to undergo an 
electromagnetic interaction and of course fulfills 

E n' I'm' _|_ br (~\ r )\ 
U nlm ' °nlm • 



em 



nlm 

n'l'm' 

The total probability for a pionic atom to suffer an electromagnetic collision per unit 
length is then given by 

n em _ 1 _ em f]Q\ 

Pnlm — \ em ~ A nlm > \ LO ) 
A nlm A 

where X^ m is the mean free path before an electromagnetic interaction takes place. 
Exploiting the completeness of the eigenstates of the Coulomb Hamiltonian the total 
electromagnetic cross sections can be calculated directly jTHJUB] and not just via (fT2|) §. 

The electromagnetic cross sections have been obtained with different approaches 
in [SlEllZj- We will devote section to discussing the different break-up probabilities 
they lead to. 

To get an insight into the magnitude of these interaction probabilities we show 
in figure Elthe average values of the annihilation, ionization, excitation and de-excitation 
probabilities per unit length. The average is taken over the even z-parity states (i.e., 
/ — m even) for fixed n. The atoms are created in even z-parity states (I — m — 0) 
and the transitions to odd ^-parity ones are strongly suppressed. The figure shows the 
probabilities using the coherent (interaction with the atom as a whole) contribution of 
the Born2 set of cross sections. This cross section set will be described in section |3J 
Any other choice among the cross section sets described in section HI would lead to very 
similar results. The averages are defined as 

= n(n + l)/2 ^ ' (14) 

lm 

P " r = n(n + 1)/2^" L ' (15) 

V ;/ lm 

lm n'<n,l'm' 

^^ s^y^ ip^rn-P^- E PnLA , (17) 
^ >' lm \ n'<n,l'm' / 

where n(n + l)/2 is the number of even z-parity states for a given n. 

§ Note, however, that this is strictly true only within the framework of the sudden approximation 
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n 



Figure 2. Annihilation, ionization, excitation and de-excitation probabilities per unit 
length, averaged over I and m quantum numbers according to equations (|14|) . Q15JI. 1)16(1 
and 1)17(1 , as a function of the principal quantum number n. 



Pionium Evolution in the Target 
To simulate the evolution of a pionic atom we use the following algorithm: 

(i) We generate a laboratory momentum P, an initial set of quantum numbers and an 
initial position R for the atom as described in subsection 12.11 

(ii) We generate a free path according to: 

p(x) dx = -}—e~ x/Xnlrn dx (18) 

where \ n i m = (l/^nim + V^n?m) _1 is the mean free path before either an 
electromagnetic interaction or the annihilation takes place. 

(iii) We displace the atom by the distance x: 

p 

R> = x— + R. (19) 

(iv) We determine whether the atom has been annihilated, excited (or de-excited) in a 
discrete collision, or broken up. The relative weights of the respective branches of 
the evolution are given by the probabilities of equations (jSJ), (fTUjl . and (fTTf . 

(v) If the atom has been scattered and suffered a discrete transition we return to step 
(ii) using the new quantum numbers n', I' and m' and the new position R 1 as the 
initial values. 

More details on this model may be found in [TTj . 
3. Break-up Probability Calculation 

In principle, the break-up probability calculation of pionium should be straightforward 
once we have established the Monte Carlo model. The rest would be a matter of 
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generating an atom sample and computing how many of them break up in the target. 
However, two main difficulties arise when trying to implement the algorithm. 

The first difficulty is due to the presence of an infinite number of atomic bound 
states in the calculations. Clearly, only a finite number of states can be taken into 
account in the simulation of the evolution of pionium. In our calculations we have 
imposed a cut on the states with n < n max . This would not pose a serious problem 
if the atoms, being created mainly in very low n states, could not get highly excited. 
Unfortunately, excitation to ever higher lying bound states constitutes a major branch 
in the evolution of pionium. As a consequence we cannot directly calculate the break-up 
probability as outlined in the previous paragraph. 

The other difficulty lies in the fact that for some of the cross section sets to be 
studied in section 0J the break-up cross sections have not been calculated. In this case 
it is imperative to find an indirect way to compute the break-up probability. 

We have discussed in the previous section that pionium terminates its evolution in 
the target by being either annihilated or broken up. However, the atom can also leave 
the target in a bound state. This would happen if one of the generated free paths in 
the Monte Carlo procedure carries it to a position outside the target. The break-up 
probability (Pbr), the annihilation probability (P an h), and the probability to leave the 
target in a discrete state (Pdsc) are related by: 



This equation allows us to compute the break-up probability indirectly. 

3.1. Computation Difficulties Due to Physical Characteristics of the Problem 

The probability to generate an atom in a specific shell decreases as 1/n 3 . This means 
that the number of atoms created with n > 4 is very small. If the atoms could not 
get excited to states with large n, we could safely solve the evolution system by setting 
n max > 4. However, as we saw in figure El the atoms have a tendency to be excited, as 
n increases, rather than being annihilated or ionized. 

Hence we expect a significant fraction of atoms excited into n > n max shells even 
for large values of n max . The probability of an atom in a nlm state to be excited into a 
state beyond the cut, i.e. with n > n max , is given by 



where we have used (jllj) and (fT3|) . However, once the atom jumps into one of 

these states we loose control over it and we have to stop its evolution. 

To analyze the change of the Monte Carlo results with n max we have modeled the 
evolution of a sample of atoms by changing n max from 7 to 9. We observed three main 
effects: 

• The fraction of annihilated atoms (P an h(^ < nmax)) does not change significantly. 




(20) 




nlm 



(21) 
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• The portion of atoms leaving the target in discrete states (Pd sc (^ < n maK )) changes 
only slightly. 

• The fraction of dissociated atoms (Phr{ n < n m3X )) changes significantly. 

This effect can be understood by checking the dependence of the annihilation, the 
discrete and the break-up probabilities on n, the principal quantum number of the 
state from which the atom was annihilated, broken-up, or in which it left the target. 
In figure El we show the result of the Monte-Carlo simulation with n max = 8 for a sample 
of one million atoms using the Born2 cross section set that also includes cross sections for 
the ionization (refer to section HJ. For the annihilated atoms we can see that P an h(^) is 
negligible for values of n > 4. The Pd sc (^) dependence also shows a fast, but less drastic, 
decrease with n. Only the solution for the states with n = n max — 1 or n = n max — 2 
is unstable under variation of n max . For n max = 8 this is a small contribution to the 
total Pdsc value. Finally, -Pbr(^) decreases very slowly with n, showing that there is a 
significant fraction of atoms broken up from states with n > ri max . The probability of 
an atom to be excited into such a state with n > n max is also shown as a function of the 
principal quantum number of the last state before the excitation. Obviously, this effect 
is non- negligible. 

For the cross section sets without break-up cross sections, we can calculate directly 
only the total probability for all electromagnetic processes and the probabilities for 
discrete transitions to states with n' < n max . In these cases we can therefore not 
distinguish whether an atom has been broken up or excited into a state with n' > n max , 
that is, we can only determine the combination of probabilities 

n'>n max _|_ br _ em _ n'l'm' (22) 

"nlm ' "nlm "nlm / j "nlm ' \ I 

n'<n max _,l'm' 

This is, of course, equivalent to (jUJ), but in this case p^ m is unknown. Thus, for these 
cross section sets not even the break-up probability for low n states could be directly 
calculated and we are forced to use the procedure described below. 

3.2. Calculation Procedure 

Based on the fast decrease of P arm and Pdsc as a function of n we can assume that 
almost every atom excited to a state n > n max will be eventually broken up. This will 
be true even though the excitation probability per unit length of a given bound state 
is significantly larger than the break-up probability per unit length. We can explain it 
as follows. The mean free path of the excited atoms strongly decreases with increasing 
n. For n ~ 8 the mean free path is < 0.1 /im. An excited atom will thus interact 
many times within a very short distance. In every scattering the atom will have some 
small probability to break up, thereby terminating its evolution. In summary, the most 
probable evolution of an atom that has been excited to any state with n > 4 is a 
sequence of excitations (and less frequent de-excitations) terminated by break-up. 
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Figure 3. Probabilities of finishing the evolution in a discrete state (a), by 
annihilation (b), or by ionization (c) as a function of the parent state's principal 
quantum number. In (d) we show the probability for an atom in an state n to be 
excited into a non-controlled state with n > 8. The results are for pionium in a 
95 //m Ni target and the lifetime is assumed to be 3 • 10 -15 s. 

However, while we are neglecting the atomic annihilation from states with n > 8 
and thus setting P ann = P an h(^ < 8), we can estimate Pd sc (^ > 8) by means of a fit to 
the Pdsc(^) histogram as recommended in jS] 

P dsc (n>8) = ^ + A (23) 
Hence, taking into account ()20j) we obtain 

-Pbr = 1 — Pise — -Panh, (24) 

where Pd sc consists of two parts, 

Pdsc = Pdsc(n < 8) + P dsc (™ > 8), (25) 

of which Pdsc(^ < 8) is computed directly and Pd sc (™ > 8) is calculated from (j2HJ). In 
this manner we can calculate the break-up probability even without ionization cross 
sections as input. 

In table ^ and in figure 0] (top left) we show a few results for the probability for 
different lifetime values in a 95 /im Ni target. The target choice coincides with that of 
the DIRAC experiment. We observe that the result of Pd sc (^ > 8) adds only a small 
correction. In figure E] we also show the ionization and annihilation distributions as a 
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Table 1. Results for the different probabilities defined in (|20|l . as calculated with 
the Born2 cross section set for a sample of ten million pionic atoms in a 95 /im thick 
Nickel target. 



t[1(T 15 s] 


-Pbr 


-Panh 


P dsc (n < 8) 


P dsc (n > 8) 


1 


0.2976 


0.6527 


0.0491 


0.0006 


2 


0.3951 


0.5287 


0.0754 


0.0008 


3 


0.4599 


0.4451 


0.0941 


0.0009 


4 


0.5062 


0.3848 


0.1080 


0.0010 


5 


0.5408 


0.3392 


0.1190 


0.0010 


6 


0.5681 


0.3029 


0.1279 


0.0011 


7 


0.5901 


0.2740 


0.1348 


0.0011 



function of the target position, and finally the creation position for those atoms that 
managed to emerge from the target in a bound state. As emphasized in subsection 13. lj 
with increasing n only the atoms very near the target end will be able to leave the target 
in a discrete state. 

4. Cross Sections Sets 

In our calculations of the break-up probability we employed three different sets of cross 
sections. The first two have been calculated in the framework of the Born approximation. 
We assign the labels Bornl to the calculations made in reference [3] and Born2 to those 
of [E]. The two sets differ in four main points: 

• The Bornl set neglects the contribution of incoherent scattering (collisions leading 
to an excitation of the target atom), thus considering the coherent interaction only 
(collisions with the target atom as a whole), i.e. the leading term. By contrast the 
Born2 set accounts for target excitations. 

• The Bornl set uses Moliere's parameterization ^H] for the Thomas-Fermi equation 
solution as the target atom form factor of the pure electric interaction, whereas the 
Born2 set takes electron orbitals determined numerically within the Hartree-Fock 
framework for the same purpose. The Thomas- Fermi- Moliere parameterization of 
the atomic form factor is accurate for low momentum exchange, but gives a small 
excess for harder scattering. 

• The Bornl set considers the sudden approximation (no recoil energy for the target 
and the pionic atom) and neglects the energy difference between the initial and the 
final state, while the Born2 set accounts for these two effects. 

• Finally, Born2 set also considers the effect of magnetic and relativistic terms. 

In principle it has been concluded jE] that accounting for second order effects 
like the magnetic terms of the Hamiltonian, the recoil energy of the atoms, or the 
relativistic terms generally leads to an overall decrease of the sudden approximation pure 
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Figure 4. Top left: The break-up, annihilation, and discrete probabilities as a 
function of lifetime. Top right: Break-up and annihilation position distributions. 
Bottom: Creation position of those atoms that leave the target in a bound state 
(and contribute to Pdsc)- Note that as n increases, only the atoms very near the target 
end can escape from it. All three plots refer to a 95 /xm Ni target. In the last two, the 
lifetime is assumed to be 3 • 10~ 15 s. 



electrostatic coherent cross section value due to destructive interference with the leading 
orders. Moreover, employing atomic orbitals obtained in the Hartree-Fock approach for 
the form factor used to compute these cross sections leads to lower values than those of 
the Bornl set since the Moliere parameterization of the solution to the Thomas- Fermi 
equation is excessive for mean and large values of the photon momentum transfer. This 
last issue leads to discrepancies that increase for large n states and decrease for large 
Z target atoms. The difference appears to be balanced by neglecting the incoherent 
contribution to the cross section in the Bornl. This results in a systematically smaller 
ground state cross section of Bornl set whereas for larger n Bornl cross sections are 
larger (up to ~ 10% discrepancy) or compatible with BornS results. We can observe the 
comparison for three target materials in figure 03 and we shall analyze the disagreement 
in the break-up probability resulting from this effect. 
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Figure 5. In the left column we compare the Bornl and Born2 cross section sets. 
The middle column shows a comparison of the Bornl and Glauber sets and finally 
on the right we compare the Born2 and Glauber sets. The comparison is made for 
Titanium (Z = 22), Nickel (Z = 28) and Platinum (Z = 78). The plots refer to total 
electromagnetic cross sections averaged over m for even z-parity states. 



Finally we have also used a set of cross sections where the Glauber formalism has 
been applied to calculate the coherent contribution to the cross section value. The details 
are given in [7j. This calculation technique accounts for multi-photon exchange in the 
pionium-target atom collision. Contrary to what one would expect the consideration 
of more than one photon being exchanged diminishes the values of the cross sections 
due to a destructive interference of the n-photon exchange contributions (this happened 
also when accounting for magnetic terms in the Born2 set). The leading order of the 
Glauber result matches the sudden approximation of the Born cross sections (since both 
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neglect the difference between the initial and the final state energies). However, this 
cross section set uses a parameterization for the target atom form factors similar to the 
ones used in Born2. This explains the disagreement with respect to the Bornl and the 
agreement with Born2 set for low Z targets, as can be seen in figure EJ The corrections 
due to multiphoton exchange are important for large Z targets [Zj and this explains the 
large discrepancies obtained for Platinum. 

5. Results and Conclusions 

After the discussions of the previous sections, we finally present the results of the break- 
up probability calculation. The DIRAC experiment has the possibility of choosing 
between several targets. The design of these different targets has been made to achieve 
the maximum break-up probability resolution in different lifetime ranges. Large Z 
targets with larger interaction cross sections are better suited for small lifetime values 
whereas lower Z materials are more sensitive to larger lifetime values. Three of these 
targets are the Pt 28 /zm target, suitable for lifetime ranges r < 1 • 10~ 15 s, the Ni 
95 /im target for r ~ 3 ■ 10 -15 s and the Ti 251 jum target for r ~ 4 ■ 10~ 15 s. The 
target thickness was chosen so as to have the same radiation length and hence equivalent 
multiple scattering effects for all three targets. The Nickel target constitutes DIRAC's 
main target with which 90% of the data have been collected, as it is optimal for the 
theoretically predicted lifetime value. 

In figure El the break-up probability curves are shown for these three targets. The 
calculation has been carried out for samples of ten million events, with a statistical error 
less than 0.08%. 

One can clearly see that for the Ti and Ni targets the Glauber and Born2 sets 
lead to similar results whereas the Bornl set shows a ~ 8% disagreement. For the 
large Z target (Pt) both Bornl and Born2 are biased toward large values. In this case, 
the multi-photon contributions to the total cross sections are not negligible. In any 
case the discrepancies between the break-up probability results are at the level of the 
discrepancies between ground state cross sections and much smaller than the differences 
between cross sections sets for medium or highly excited states. We can understand this 
based on the fact that the probability for the atoms to leave the target in a discrete 
state other than the ground state and maybe the first excited shell Pdsc(^ > 2) is of 
the order of, or smaller than, 5%. Hence, even large uncertainties in this magnitude 
(up to 10 — 15%) lead to very small changes in the break-up probability result. Only 
discrepancies in the ground state population and maybe the first excited shell, where 
most of the atoms are created, would lead to significant differences between the break-up 
probability results of the different sets. 

Graphically we can view the atom as a balloon being inflated in every collision with 
the target. The different sets will lead to similar results of the size increase rate as long 
as the atom remains in a low excited state. However, as the atom grows (inflates) it 
will no longer be able to advance as easily in the target due to its large size and will 
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Table 2. Comparison of the break-up probability results for the Ti 251 //m, the Ni 
95 fim and the Pt 28 fim targets. The lifetime value is assumed to be 3 • 10 -15 s in the 
calculations. 



Target 


pBornl 
Mr 


pBorn2 
Mr 


p Glauber 
Mr 


pBorn2 
| Mr 

pBornl 
Mr 


p Glauber 
| Mr 

pBornl 
Mr 


p Glauber 
■ Mr 

pBorn2 
Mr 


Ti 


0.3026 


0.3249 


0.3232 


-7.4% 


-6.8% 


0.5% 


Ni 


0.4425 


0.4599 


0.4555 


-3.9% 


-2.9% 


1.0% 


Pt 


0.7137 


0.7196 


0.6924 


-0.8% 


3.0% 


3.8% 



finally break-up (explode). Large discrepancies in the excitation and break-up rate of 
the excited atom will not be important given that the mean free paths for the excited 
states are very small compared to the target dimensions. 

In summary, we recall that the high precision measurement attempted by the 
DIRAC collaboration also requires an accuracy to better than 1% in our theoretical 
break-up probability calculations. We note that the seemingly large discrepancies 
among our different cross section sets particularly for pionium transitions starting from 
highly excited states do not lead to significant differences in the theoretical break- 
up probabilities. The discrepancies between break-up probabilities are coming almost 
entirely from differences in the cross sections for the lowest lying states, where both 
the atomic structure of the target and the multi-photon transitions need to be treated 
as accurately as possible. This challenge, however, has already been mastered in our 
previous works [BJ EE] where we showed that the required 1% accuracy can be achieved 
with our calculations, albeit only with the Born2 and the Glauber sets for low Z and 
with the Glauber set for large Z targets. The important conclusion of the present 
investigation is the finding that the (infinitely many!) highly excited states of pionium 
do not limit the validity of our approach even though we can include explicitly only a 
moderate number of these states in our simulations. 

Acknowledgments 

We would like to thank B. Adeva, L.L. Nemenov, L. Tauscher and D. Trautmann for 
their support. We are also indebted to M. P16, J.J. Saborido and A.V. Tarasov for their 
invaluable help. 

References 

[1] Colangelo G, Gasser J and Leutwyler H 2000 Phys. Lett. B 488 261. 

[2] Adeva B et al 1995 Lifetime measurement of tt + tt~ atoms to test low energy QCD predictions 

CERN/SPSLC 95-1 (Geneva: CERN); http://www.cern.cri/DIRAC 
[3] Nemenov L L 1985 Sov. J. Nucl. Phys. 41 629. 

[4] Gortchakov O E, Kuptsov A V, Nemenov L L and Riabkov D Yu 1996 Phys. of At. Nucl. 59 2015. 

[5] Afanasyev L G and Tarasov A V 1996 Phys. of At. Nucl. 59 2130. 

[6] Halabuka Z, Heim T A, Trautmann D and Baur G 1999 Nucl. Phys. 554 86. 



A Monte Carlo Calculation of the Pionium Break-up Probability 



16 





























« 










.•B A 














W A 












fA 














Of* 

/A 








Lriauoer 














Born 2 




A 


251 


urn Ti 




















Born 1 




- , , , , 




, , , , 


, , , , 


, , , , 


, , , , 




, , , 



x[10" 15 s] 

























































' 


Glauber 
Born 2 






95 p 




a 




Rom 1 




" , , , , 




, , , 




, , , 


, , , 




, , , 



x[10" 15 s] 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 










































































Glauber 












• 








28 (J 


LmP 


t. 




Born 2 














Born 1 




- , , , , 




, , , 




, , , 


, , , 




, , , 



x[10" 15 s] 



Figure 6. The break-up probability results for the three cross section sets and the 
three target materials. 
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